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ABSTRACT 

We use the “Dark Energy and Massive Neutrino Universe” (DEMNUni) simulations 
to compare the constraining power of “sufficient statistics” with the standard matter 
power spectrum on the sum of neutrino masses, M v = In general, the power 

spectrum, even supplemented with higher moments of the distribution, captures only 
a fraction of the available cosmological information due to correlations between the 
Fourier modes. In contrast, the non-linear transform of sufficient statistics, approx¬ 
imated by a logarithmic mapping A = ln(l + 6), was designed to capture all the 
available cosmological information contained in the matter clustering; in this sense it 
is an optimal observab le. Our analysis takes advantage of the recent analytical model 
developed by Carron et al. 2014) to estimate both the matter power spectrum and 
the A-power spectrum covariance matrices. Using a Fisher information approach, we 
find that using sufficient statistics increases up to 8 times the available information on 
the total neutrino mass at z = 0, thus tightening the constraints by almost a factor of 
3 compared to the matter power spectrum. 
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methods, cosmology : cosmological parameters, methods: simulations 


1 INTRODUCTION 

The interest for neutrino science has been stimulated over 
the last decade by solar, atmospheric and accelerator ex¬ 
periments that showed, through the observations of the so- 
called “neutrino flavour oscillations” that neutrinos have a 


finite mass ( 

Ahmed et al. 2004 Eguchi et al. 2003 Araki 

|et al. 2005 

McKeown & Vogel 2004). These experiments 

are, however, only sensitive to mass square differences be- 


tween neutrino mass eigenstates, leaving the knowledge of 
the total neutrino mass as one of the great unsolved prob¬ 
lems of modern physics. Cosmological probes, such as the 
cosmic microwave background (CMB) or large scale struc¬ 
tures (LSS) in the Universe have the highest experimental 
sensitivity to date to the absolute neutrino mass, as mas¬ 


sive neutrinos influence structure formation (Dolgov 2002 
Lesgourgues fe Pastor||2006 20121. 


The epoch at which massive neutrinos became non- 
relativistic changes the time of matter-radiation equality, in¬ 
creasing slightly the size of the sound horizon at recombina¬ 
tion and thus changing the position of the CMB anisotropy 


peaks. In combination with Baryon Acoustic Oscillations 
(BAO), measurements from the Planck satellite have pro¬ 
vided the tightest constraints on the total neutrino mass 
M„ < 0.17 eV at 95% CL (Planck Collaboration et al. 20151. 


However this result is sensitive to the value of the Hubble 
parameter and including other external measurements leads 
to the weaker, more conservative limit of M v < 0.23 eV. 
Neutrinos also impact the growth of structure by suppress¬ 
ing the small-scale matter density fluctuations. Thus mea¬ 
surements of the total matter power spectrum can improve 
constraints on the neutrino mass in a complementary way. 
The Sloan Digital Sky Survey (SDSS) DR8 LRG angular 
power spectrum with the WMAP7 data and a HST prior 
on the Hubble parameter gives M v < 0.26 eV (95% CL) 
(de Putter et al. |2012[ ) . More recently |Beutler et al.| ( |2014[ ) , 
using the Baryon Oscillation Spectroscopic Survey (BOSS) 
CMASS Data Release 11 combined with measurements from 
the Planck satellite (without the A^-lensing signal), weak 
lensing data and BAO constraints, found M„ = 0.36 ± 0.10 
eV (68% CL). 
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These bounds mainly use structure formation data up 
to the translinear regime ( k max ~ 0.15 /iMpc^ 1 ), but 
even on such large scales some non-linear contamination 
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is present and must be taken into account. Moreover, fu¬ 
ture galaxy redshift surveys such as the Dark Energy Sur- 
ve >0(DES), the Large Synoptic Survey Telescop^] (LSSTL 
Euclitj^] or the Wide Field Infra-Red Survey Telescopcj^] 
(WFIRST), are expected to constrain the matter power 
spectrum at percent level precision at least on scales of 
k ~ 0.1 — 1 hMpc -1 ( Hearin et al.|2012 i. 

Yet, the power spectrum, as well as higher order statis¬ 
tics, lose their effectiveness beyond linear scales due to the 
emergence of coupling among Fourier modes induced by the 
non-linear gravitational evolution (Rimes & Hamilton 2005 


Neyrinck et al. 20061. It means that even if these future 


surveys can probe more deeply the non-linear regime, the 
amount of accessible information on the neutrino mass, or 
in general on other cosmological parameters, will be drasti¬ 
cally reduced. Non-linear transformations, such as the log¬ 
arithmic mapping ( Neyrinck et al.|2009 l or variants thereof 
(Seo et al. ||201 1 Joachimi & Taylor 20111 were introduced 
specifically to recapture this hidden information. |Carron 


Szapudi (2013), using perturbation theory, shows that this 


logarithmic transformation, A = ln(l+<5), approximates well 
the exact “sufficient statistics” for a continuous field, mean¬ 
ing that it captures, by design, all the available cosmological 
information from the matter field. With the data analysis of 
future cosmological surveys in mind, investigating the con¬ 
straining power of this new observable via IV-body simula¬ 
tions with a massive neutrino component is necessary. Our 
goal is to extract all possible information on cosmological 
parameters, especially, the neutrino mass, the focus of the 
present work. 

We proceed as follows. Section[2]presents the DEMNUni 
simulations. The analytical model used for the estimation 
of the matter power spectrum and the A-power spectrum 
covariance matrices is described in Section [3] Making use of 
this model, we forecast the Fisher information on the total 
neutrino mass for both the non-linear transform A and the 
matter power spectra. We summarize and conclude with a 
discussion in Section [4] 


2 SIMULATIONS 


The DEMNUni simulations, presented in [Carbone et al.| 
|2015[ ); [Castorina et al.| |2015|, are the largest N-body sim¬ 
ulations to date with a massive neutrino component treated 
as an additional particle type. At present, they are charac¬ 
terised by a baseline ACDM cosmology to which neutrinos 
are added with a degenerate mass spectrum but different 
total masses. In the near future, the DEMNUni set will in¬ 
clude also an evolving dark energy background, with differ¬ 
ent equations of state w , in order to study the degeneracy 
between M v and w, at the non linear level. 

The DEMNUni simulations include only Cold Dark 
Matter (CDM) and neutrino particles, and have been 
performed using the tree particle mesh-smoothed parti¬ 
cle hydrodynamics (TreePM-SPH) code gadget-3 |Springel| 
(20051, specifically modified in Viel et al. (20101 to account 


1 http://www.darkenergysurvey.org/ 

2 http://www.lsst.org/lsst/ 

3 http://sci.esa.int/euclid/ 

4 http://wfirst.gsfc.nasa.gov/ 


for the presence of massive neutrinos. These simulations 
have been run on the Fermi IBM BG/Q supercomputer 
at CINECA, ItaljQ They contain (2048) 1 2 3 4 CDM particles 
and (2048) 3 neutrino particles, in a comoving cube of vol¬ 
ume V = 8h _3 Gpc 3 , and are characterised by a softening 
length e = 20 Kpc /h. For each simulation, 62 outputs have 
been produced, logarithmically equispaced in the scale fac¬ 
tor a = 1/(1 + 2 ), in the redshift interval 2 = 0 — 99, 49 of 
which lay between z — 0 and 2 = 10. 

The DEMNUni set is composed by four cosmologi¬ 
cal simulations sharing the same baseline cosmology con¬ 
sistent with the cosmological parameters estimated by the 
first Planck data release (Planck Collab oration et al.|2014 |: 
fl m = 0.32, 12a = 0.68, Ho = 67kms _1 Mpc“ , Tib = 
0.05, andn s = 0.96. The only difference between the four 
simulations is represented by the total neutrino mass M„. 
One of the simulations is characterised by a pure ACDM 
cosmology without neutrinos, and is used as a reference. 
The other three simulations correspond to a ACDM cos¬ 
mology plus three degenerate massive neutrinos with M„ = 
0.17, 0.3, 0.53 eV. As massive neutrinos modify the shape of 
the power spectrum during the cosmic evolution, a special 
care has been taken in setting the initial power spectrum 
in each simulation. The power spectrum used to set up the 
initial conditions has been obtained from CAMI^] adopt¬ 
ing the same primordial scalar amplitude for all the four 
simulations. On one hand, this ensures that the normaliza¬ 
tion of the power spectrum in the four runs is consistent 
with Planck, on the other hand it guarantees that the shape 
of the initial power spectrum is consistent with the consid¬ 
ered neutrino mass. Given the assumed Planck cosmology, 
the particle number, and the volume of the simulations, the 
CDM mass resolution is about 8 x 10 10 /i _1 Mg, adjusted ac¬ 
cordingly to the value of the neutrino particle mass, in order 
to hold n m = 0.32 fixed. 


3 RESULTS 

3.1 Power spectra and covariance matrices 

The power spectra are computed using the standard estima¬ 
tor defined as 

P(fc) = ^^ \6(k')\ 2 , (1) 

k' 

where V is the survey volume and the sum runs over the 
Nk Fourier modes associated to the fc-th power spectrum 
bin. We measure the real-space power spectra of both the 
density, P, and of the non-linear transform A, Pa, defined 
as 

A = ln(l + 5), 8 = THlI ( 2 ) 

P 

on a 512 3 grid with the nearest-grid point density assign¬ 
ment deconvolving for the pixel window over the range 
0.0035 < k < 0.4 hMpc -1 . The power spectrum covariance 
matrix is defined as 

Cov»j = <£(&)£(**)> - (P(fci)>(P(fcj)>. (3) 


5 http://www.hpc.cineca.it/ 

6 http://camb.info/ 
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Table 1. Contributions to the covariance matrix from “super¬ 
survey” and “intra-survey” modes respectively as a function of 
redshift. The amplitude of both parameters are estimated using 
the above definitions of and rf s - 


Table 2. Information gain on neutrino mass for different red- 
shifts. It represents the ratio between the elements of the power 
spectrum and A-power spectrum Fisher matrix. This gain can be 
interpreted as an effective gain in the survey volume. 


Redshift 

°SS 

a lS 

2 = 0 

2.3 x 10~ 6 

1.4 x 10~ 6 

2 = 0.5 

1.4 x 10 -6 

6.0 x 10 -7 

2 = 1 

9.0 x 10 -7 

3.2 x 10 -7 

2 = 1.5 

5.7 x 10~ 7 

2.0 x 10- 7 

2 = 2 

4.1 x 10" 7 

1.3 x 10 -7 


We follow Carron et al. (2014) who developed an approxi¬ 


mate form of the latter in the mildly non-linear regime based 
on previous studies from 7V-body simulations ( Neyrinck] 
2011 Mohammed &; Seljak|2014 ) 


COV ij — Sij 


2 P(ki) 2 
N ki 


+ % P{ki)P{kj 


(4) 


where the first term corresponds to the Gaussian covari¬ 
ance and the second term approximates the shell-averaged 
trispectrum of the held. Carron et al. (2014) have shown 
that the parameter cr 2 lin can be interpreted as the min¬ 
imum variance achievable on an amplitude-like parame¬ 
ter. It can be further decomposed into two contributions 
a rnin = &ss + a is> where the first term is due to the corre¬ 
lation between large wavelength “super-survey” modes with 
small scales and the second term corresponds to coupling 
between small scales or “intra-survey” modes. 

As we consider the local density fluctuations, S = 
defined with respect to the local observed density (i.e 
the global density fluctuations divided by the background 
mode), we have an additional contribution to the covari¬ 
ance matrix that can be entirely absorbed within cr| s . The 
latter can thus be expressed as agg = (26/21) 2 oy, where 
following Takada & Hu (2013) 


<jy = / P L (k)W 2 (k) 


drk 

(27r) 3 


(5) 


In this equation, W (k) is the Fourier transform of the survey 
volume window function and P L is the linear power spec¬ 
trum ( |Carron et al.|[2014| |Takada fe Hu||2013| for details). 
Here W[x ) is defined as a spherical top-hat on the volume 
V. 


The contribution to the covariance matrix coming from 
ajg could be more complex, as it was only tested at 2 = 0 
( Mohammed fe Seljak|2014 l. However, Carron et al. (20141 
derived an analytical approximation based on the hierarchi¬ 
cal Ansatz ( Peebles||1980 Fry||1984 Bernardeau||1996 1 and 
found that 


2 

a is 


p{k ;; ax) (4 Ra + 4 Rb ). 


( 6 ) 


We assume this expression holds for all redshifts and use it 
to estimate the a 2 s contribution. Our results are presented 
in Table [I] and from it a 2 lin is calculated using the above 
definition. To quantify its Fisher information content, we 
need an estimation of the A-covariance matrix as well. We 
choose to adopt a Gaussian description for the A field 


Cowt 


2 

Nk 


PA{ki)PA(kj)5ij. 


(7) 
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Figure 1. Logarithmic derivatives with respect to at 2 = 0 
of the power spectrum (black lines) and the A-power spectrum 
(red lines) at two different fiducial neutrino masses. 


Neyrinck (2011) have shown this form is valid to a very 
good approximation for k < 0.4 /iMpc” 1 , scale at which 
the covariance matrix of Pa starts to have non-negligible 
off-diagonal elements. We thus make our measurements up 
to this resolution limit, letting us safely assume that Pa 
has a diagonal covariance matrix over the full range of k 
considered here. 


3.2 Fisher information 


Given a set of parameters a,/3 ,..., the Fisher matrix of the 
matter power spectrum is defined as 


F aP = £ 

ki,kj <k n 


dP(ki) idPfe) 

~fcT CoVij ~W~- 


( 8 ) 


of the power spectrum can be written as 


Carron et al. (J2014) have shown that the information content 

(9) 


FaB = F n g ~ a. 


rpG rpG 
2 r aln Ao^ln 


where 


pG _ 

-F rv R 


/ 


1 + °Ln(S/N) a a’ 

d\nP(k) d In P(k) 


da 


d/3 


-w(k) dink, 


( 10 ) 


with w(k) = V k 3 /(2n) 2 . The discrete sums have been re¬ 
placed by integrals using the fact that the number of modes 
N k is approximately the surface of the shell used for the 
bin averaging divided by the distance element between two 
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discrete modes AT*, — V (4tt k 2 dk) /(2n) 3 . Moreover, in Equa¬ 
tion [ 9 ] by analogy to Carron et al. (20141, we have intro¬ 
duced, for notational convenience, a nonlinear amplitude 
parameter, InAo, defined such as d\ n A a P(k ) = P{k). This 
parameter corresponds to the initial amplitude a£ in the 
linear regime and at 2 = 0. We further define the Gaussian 
signal to noise as 

(S/N) 2 g = J w(k) dink . (11) 

Using Equation [7] the Fisher information in the A-power 
spectrum is given by = F„p A , where 


F?i A = 


/ 


d In Pa (k) d In Pa (k) 


da 


d/3 


w(k) dink. 


( 12 ) 


In this analysis we focus on the information about the neu¬ 
trino mass, we thus consider a = /3 = In (M v ) and hold fixed 
all the other cosmological parameters. 

The derivatives in Equations [TO] and [12] are calculated 
numerically by comparing the power spectra among sim¬ 
ulations with different total neutrino masses. The deriva¬ 
tives at 2 = 0 for both the density and the A fields are 
shown in Figure |7] at the minimum and maximum neutrino 
masses considered: M v = (0.17 + 0.30)/2 = 0.235 eV and 
M v = (0.30 + 0.53)/2 = 0.415 eV. 

The results are presented in Figure [2] which shows the 
cumulative Fisher information on the logarithm of the neu¬ 
trino mass using both the power spectrum (black lines) and 
the A-power spectrum (red lines) at two different redshifts: 
2 = 0 (left panel) and 2 = 1 (right panel). The Gaussian 
part of the matter power spectrum is given for comparison 
by the blue solid line in the case M v = 0.415 eV. At our 
resolution limit, using sufficient statistics increases by up to 
~ 8 the available information, leading to error bars that are 
almost a factor of 3 smaller compared to the matter power 
spectrum at 2 = 0. The evolution of the information gain 
with the redshift is presented in Table [2] As was shown pre¬ 
viously by Wolk et al. (20141 in the case of projected field, 


the gain is larger at small redshift where the non-linearities 
are stronger. The fiducial value of the neutrino mass used 
to estimate the derivatives has almost no impact on the fi¬ 
nal gain with, however, a small tendency to have lower gain 
at higher neutrino mass as it can be seen on Figure [2] or 
in Table [2] Our results show that even at 2 = 2, the suf¬ 
ficient statistics unveils up to ~ 2 times more information 
compared to the matter power spectrum. These conclusions 
directly depend on the value of and a change of 10% 

of the latter will result in a change of about ~ 5% in the 
information gain. 

It can be seen on Figure [5] that in the regime 0.025 < 
fc < 0.1 AMpc” 1 , Pa is less optimal than the matter power 

spectrum. This is due to the bias between Pa and P that 

2 

could be, in a first approximation, written as Pa = e A Pu„ 


where a\ is the variance of the A-field (Neyrinck et al. 


20091. As this bias depends on the cosmology, it shifts the 


0-crossing scale of the derivatives of Equation [12] from a 
scale k ~ k m i n for the matter power spectrum to k ~ 0.03 
/iMpc -1 (see Figure [l]), resulting in a knee in the informa¬ 
tion. Moving to larger fc, the derivatives increase again in 
absolute value and Pa starts again to perform better than 
the matter power spectrum. 


4 DISCUSSION 

It has been extensively shown that the power spectrum is 
not the optimal observable to constrain cosmological param¬ 
eters in the mildly non-linear regime as the information it 
contains saturates at a finite plateau thus contradicting the 
naive Gaussian expectation (Rimes fc Hamilton|2005] 


2006 


Neyrinck et al. |[2006| [Neyrinck fc Szapudi||2007[ |de Putter| 


et al. 


2012 


Takada & Hu 20131. Sufficient statistics have 


been introduced to overcome this issue and recapture, by 
design, in their power spectrum most of the available infor¬ 
mation. We have demonstrated, using the DEMNUni simu¬ 
lations, which correspond to a ACDM cosmology with three 
degenerate massive neutrinos with sum of masses respec¬ 
tively 0.17, 0.3 and 0.53 eV, that the power spectrum of 
the non-linear transform A outperforms by a factor ~ 8 the 
usual power spectrum, when extracting cosmological infor¬ 
mation on the neutrino mass at 2 = 0, and by a factor of 
~ 2 even at 2 = 2. This gain can be seen as an effective gain 
in the survey volume. 

Observational issues were not considered in this work, 
and refinements need to be made in order to apply the non¬ 
linear transform A to actual galaxy survey data. The first 
one is the discreteness of galaxy surveys that must be taken 
into account. The effect of the shot-noise on F a: p will be to 
change in(fc) into 


w(k) = V 


k 3 


(2^) 2 


w(k) = V 


fc 3 


( Pik) 

(2tt) 2 VP(fc) + 


\ 2 

;) , (13) 


and P(k) —> P g = b 2 P(k), where b is the linear galaxy bias. 
In the case of the “sufficient statistics” a new observable A *, 
optimized for both the non-linearities and the observational 
noise, has to be used to recapture the cosmological infor¬ 


mation (Carron & Szapudi 2014 Wolk et al. 2014). Wolk 


et al. (20151 have derived in the case of projected fields, the 
information content of A* as a function of both the infor¬ 
mation in the galaxy power spectrum and the cosmologi¬ 
cal dependencies of the bias between P and Pa- Assuming 
this approach holds in the 3-dimensional case, we consider 
the Sloan Digital Sky Survey (SDSS) LRGs sample number 
density, h ~ 5.10~ 4 fc 3 Mpc -3 with a galaxy bias of b = 2 
( Percival et al.|2007 de Putter et al.|2012 1 . At our k max , we 
fold la error bars of ~ 0.01 and ~ 0.005 using P and Pa* 
respectively leading to an information gain on M v of ~ 4 at 
2 = 0.3 (corresponding to the peak of the distribution). As 
the clustering information content is well understood in 2D, 
this simple estimation is expected to give a realistic order 
of magnitude, however further investigations about the re¬ 
lationship between the matter and the non-linear transform 
A power spectra need to be made in 3-dimensions. 

In this work, we have considered the case of the den¬ 
sity fluctuation S defined with respect to the local observed 
density as it is the case for galaxy surveys. However, in the 
case of fluctuations defined with respect to the global den¬ 
sity, for example in weak lensing surveys, our predictions 
for the information change. The contribution from the intra¬ 
survey modes becomes negligible compared to the one from 
the super-survey modes, as the latter can now be written as 
cr| s = (68/21 ) 2 ay. This gives a value of <j| s = 1.6 x 10~ 5 
and 3.0 x 10” 6 at redshifts 2 = 0 and 2 = 2, respectively. 
The expected gain in information about the neutrino mass 
using Pa becomes about factors of ~ 25 and ~ 8, respec- 
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Figure 2. Cumulative Fisher information on the logarithm of the neutrino mass using both the power spectrum (black lines) and the 
A-power spectrum (red lines) for two different fiducial values of the neutrino mass. The covariance matrices on the two observables are 
given by Equation [I] and Equation [7] respectively. The Gaussian part of the matter power spectrum in the case A& = 0.415 eV is given 
for comparison by the blue solid line. The left panel shows that using the “sufficient statistics” at z = 0 increases by ~ 8 and ~ 6.5 the 
available information for M v = 0.235 eV and M u = 0.415 eV respectively. The right panel shows the same cumulative Fisher information 
at z = 1. At this higher redshift, the gain using the “sufficient statistics” is about ~ 4 and ~ 3.5 for M„ = 0.235 eV and M v = 0.415 eV. 
The evolution of the information gain with redshift is summarized in Table [2] 


tively, compared to the power spectrum. This behavior was 
expected as the information in the matter power spectrum 
roughly saturates at l/<r)k in . In the global case, the “super¬ 
survey” modes dominate resulting in an information plateau 
which is lower compared to the local case and thus allowing 
the non-linear transform A to perform even better. 

The effect of the galaxy bias was already taken into 
account within the framework of the “halo model” for pro¬ 
jected fields ( |Wolk et al. |2014| ), as well as the bias between 
Pa * and P g . This needs to be extended to the 3-dimensional 
case. For future galaxy survey applications, it is also crucial 
to investigate the effect of redshift distortions which is left 
for a future work. 

Despite the need of a more sophisticated analysis for 
a direct application to real galaxy surveys, this work has 
demonstrated that “sufficient statistics” are expected to be 
a powerful method to improve the future constraints on dif¬ 
ferent cosmological parameters and, in particular, on the 
neutrino mass. 
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